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Abstract 

We present a calculation of the damping of an ultra-energetic cosmic 
neutrino (UHEz^) travelling through the thermal gas of relic neutrinos, us- 
ing the formalism of finite-temperature field theory. Prom the self-energy 
diagram due to Z exchange, we obtain the annihilation cross section for 
an UHEz^ interacting with an antineutrino from the background. This 
method allows us to derive the full expression for the UHEzy transmission 
probability, taking into account the momentum of relic neutrinos. We dis- 
cuss the effect of thermal motion on the shape of the absorption dips for 
different UHEz^ fluxes as well as in the context of relic neutrino cluster- 
ing. We find that for ratios of the neutrino mass to the relic background 
temperature 10^ or smaller, the thermal broadening of the absorption lines 
could significantly affect the determination of the neutrino mass and of the 
characteristics of the population of UHE;/ sources. 



1 Introduction 

One of the predictions of Big Bang Cosmology is that the Universe is filled with 
a background of neutrinos, analogous to the cosmic microwave background, but 
with a lower temperature T^q ^ 1.95 K (1.69 x 10~^ eV) and a number density 
UyQ 56 cm""^ per species E] . The direct detection of this cosmological relic 
background is extremely difficult because of the very small interaction cross- 
section of low energy neutrinos. 
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It is therefore interesting to explore the possibihty of probing the cosmic 
neutrino background (Cz/B) with ultra high energy neutrinos (UHEi/). At high 
energies, the Z resonance in the s channel for the process vv ^ X enhances 
the probability for the interaction of an UHEi/ with the Ci/B 13 El El Q IE!)- 
This process has also been proposed as a possible mechanism for generating UHE 
cosmic rays through the hadronic decay of the Z boson ( "Z-burst" mechanism j^l 
HnHHl 1121 ESI)- In that context, it has been pointed out that some features of the 
Ci/B (in particular the relic neutrino masses, energy spectra and densities) might 
be indirectly inferred from the observed spectrum of UHE cosmic rays ITK] . 

On the other hand, the resonant production of Z through vv annihilation also 
results in absorption dips in the UHEi/ spectrum. If these absorption lines can 
be observed at Earth with the appropriated resolution, their study would lead us 
more directly to the same goal: to perform relic neutrino spectroscopy, thereby 
providing us with evidence for the existence of the Ci/B and with an independent 
way of determining the absolute neutrino mass (for recent discussions on the 
subject, see [IHIIIZI)- 

Most of the work in the literature assumes that the relic neutrinos are at rest. 
For small neutrino masses, though, the average momentum can be comparable 
to the neutrino mass. In this paper, we compute the dominant contribution to 
the interaction of an UHEz/ with the Cz/B using finite temperature field theory 
(FTFT). This formalism allows us to take effects due to the neutrino background 
into account in a systematic and elegant way. 

In section |2l we evaluate the damping of a UHEz/ travelling through the Cz/B 
in terms of the imaginary part of the neutrino self energy. From the damping we 
determine the absorption probability for an UHEz/ emitted at a given redshift. 

We then present in section El some illustrations of our calculation in realistic 
physical contexts. First, the shape and position of the absorption lines in the 
spectrum of UHEz/ from interactions with the Cz/B depend on the mass of the 
neutrinos and on the type and distribution of sources for UHEz/. We explore 
various combinations of parameters to investigate the differences between the fi- 
nite temperature calculation and previous approximations and to determine the 
regimes in which those approximations break down and thermal effects become 
significant. Then we further illustrate our results in the context of relic neu- 
trino clustering for different hypothesis on the density and scale of the clusters. 
Conclusions are drawn in section jH 
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2 Damping rate and transmission probability of 
an UHEz/ 



2.1 Self-energy in the relic neutrinos thermal background 

The dispersion relation of a particle that propagates through a medium is deter- 
mined from the linear part of the effective field equation. In momentum space, 
for a neutrino with four- momentum k^^ and mass m^, it takes the form 

(^-m, -S) ^^ = 0, (1) 

where S corresponds to the retarded self-energy and embodies the background 
effects. For Dirac neutrinos, the chiral nature of neutrino interactions implies 
that, to one-loop order, 

^ = {a^ + bi)L, (2) 

where L = (1 — 75)/2 and is the velocity four- vector of the medium; in its 
own rest frame = (1, 0) and = {S, K). The coefficients a and h are complex 
functions of the scalars 

S = k.U, K = _ (^3) 

with K = \K\. In the present context S corresponds to the Feynman diagram 
of fig. d where the loop contains a relic (anti)neutrino from the thermal bath, 
with four- momentum = [Epjp), and a Z boson with a blob indicating that we 
consider its decay width to all possible channels. 

A consequence of the presence of a self-energy term in the equation of motion 
is to modify the dispersion relation of the incoming neutrino to 

Sk = Sr- i^, (4) 







u{k) 







Figure 1: Feynman diagram for the one-loop self-energy of an UHE neutrino due to a Z-boson 
exchange with an (anti-) neutrino from the relic background; the blob on the Z propagator 
indicates that we use the dressed propagator and the cut is to select the imaginary part of the 
diagram. 
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where 8k ■, Sr and 7 are functions of K. The real part, Sr, is in general not equal 
to yK^^+m? but has some additional correction reflecting the dispersive inter- 
actions that can take place in the medium, while the imaginary part corresponds 
to the damping factor, or else said, to the total reaction rate For a constant 
damping, the survival probability of an UHE neutrino travelling through the relic 
neutrino background is 

Piir) = e-^^ (5) 

function of the propagation time r. 
The damping factor is directly related to the imaginary part of the self-energy 
Sj. In the real-time formalism of the FTFT, both the propagators and the self- 
energies become 2x2 matrices. As shown in ^^1, can be expressed in terms 
of the off-diagonal elements of the self-energy matrix as: 

= I (S12 - S21) . (6) 

The vertices of the theory are doubled compared to the vacuum case; the 
subscript 1 denotes the normal vertices of the Standard Model while the ver- 
tices labelled 2 get an extra minus sign. In our case, the expressions for S12 as 
calculated from the Feynman diagram of fig. ^ is 

- ^ = 1 1% ^ Dl^{q) 7, L z 7. L, (7) 

2cos6'w J (27r)'^ 

with D2i{q) and Suip) denoting the corresponding element of the Z boson and 
neutrino propagator, respectively. One obtains S21 by interchanging the indices 
1 and 2 everywhere in eq. ((Tj). 

Since the temperature of the medium is small compared to the boson mass, 
we may discard the thermal contribution to the Z propagator. Consequently, we 
have 

D>^^{q) = 2i ImD^'^(g) e{-q.u), (8) 
= 2t \mD^^{q) Oiq.u), (9) 

where 6 is the step function and D^^^{q) is the vacuum propagator for the Z boson. 
For D'^^[q) we adopt the usual prescription around the resonance in the unitary 
norm [2011^11221: 

D'-i^) = 2 7,2 ^\r - (10) 



Here M is the Z mass and F is the total width for Z decaying to fermion pairs. 
At lowest order F can be expressed, neglecting fermion masses, as: 



V2Gf 
Gtt 



J2 ((4)' - sin^ 9y, + 2Q) sin^ , (H) 
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where Gp is the Fermi constant, N'^ = 1(3) for leptons (quarks), and Ij and Qf 
are the fermion isospin and charge, respectively [2D1- The numerical values for F 
and M are taken from ^j. We assume that the Z boson decays in vacuum. 
The propagators for the relic neutrino entering eq. are 

Suip) = 2m5ip^-ml) [7]F-ei-p.u)] i^ + m,), (12) 

S2iip) = 2711 dip" -ml) [r^^-9{p.u)] (;^ + m,), (13) 

where 

= e{p.u)UP) + e{-'p.u)u{P) (14) 

Q{P-u) 9{-p.u) 
e/3P-a + i + e^^+" + l' ^ ' 

with P denoting the magnitude of the three-momentum in the rest-frame of the 
background. The distribution of background relic (anti)neutrinos fu{P) {fu{P)) 
assumes a relativistic Fermi-Dirac form, with (3 = l/T^, and a = fi/T^y, fi being 
the chemical potential for neutrinos and T,^ the temperature of the relic neutrino 
bath. It is worth remarking here that, although relic neutrinos are not relativistic 
anymore at present time, their distribution maintains a relativistic form since 
their interactions froze out at their decoupling time, corresponding to T^d ~ 1 
MeV H. 

After replacing the expressions above in eq. expanding the delta function 
and neglecting contributions of order m'l/M'^ coming from the term proportional 
to q^qu in (fTUj) . one gets that 

+ lm[Di"'{p~k)\\e[£-E^)(l-U(P)) + 0(E^-S)U(P)]}, (16) 

where the integral is evaluated in the rest frame of the medium and we replaced 
by Ep = (P^ + m^)-*^/^ everywhere. 

The first term corresponds to the resonant production of a Z boson through 
u — u annihilation; the factor of fp{P) reflects the Pauli blocking acting on the 
antineutrinos in the background. The second and third terms correspond to the 
emission of a Z-boson respectively by the incoming neutrino {Sr > Ep) or by the 
background antineutrino {Ep > Sr). Both processes are kinematically forbidden 
and we drop them from this point on. 

2.2 Ultrarelativistic approximation 

Let us assume that the neutrino travelling through the relic neutrino thermal 
bath is ultrarelativistic and that we can neglect the background effects on its 
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energy, i.e., 8r — K. In that case, we can use the following expression for the 
damping: 

7^-2 Im6(i^,ir) = --iTr(^E,)|^^=^. (17) 

This expression was derived in jTH| for the case of a massless fermion, but it can 
be shown that it remains valid for a massive one in the relativistic limit. The term 
proportional to qaqp in the Z propagator gives contributions of order ml/M'^; if 
we discard them, then the damping rate corresponding to the uu annihilation 
process reads 

T f d'P MP) {k + p)\k.p) 

l{K) -- 



cos^^wMJ {2n)^2KEp {1 + ^){k + - 2M^{k + + M^' 

where ^ = T'^/M'^ <^ 1 and p^ = Ep. The previous formula can be conveniently 
rewritten in terms of the momentum integration of the cross-section for the pro- 
cess vv — >■ Z, weighted by the corresponding statistical factor: 

i-co p 

7(^) = y^ ^P' MP)aMP,K), (18) 



where 

a-iPK) = ^IK^ r s - (19) 

^.AP,^) ^2K^PEpJ,_ ■ (s-M^Y + ^s^- ^ 

Neglecting the chemical potential, the distribution function of the relic anti- 
neutrinos is: 

MP) = MP) = ^7v5^' (20) 

The integration variable s = {k + pY corresponds to the total energy in the 
center-of-mass and 

s± = 2ml + 2K{Ep±P). (21) 
The integral in eq. (|19|) can be done in a closed form. For irii, -C M, K, we 

get 

2V2G^TM \ 1 
a^p(P, K) = — — — <^ — — + 



X In 



2KEp il + e 4:KP{l + 0^ 

;i + 04^^(^P + P)^ - 4:M^K{Ep + P) + M^ 
1 + ^)AK^{Ep - Pf - AM^K{Ep - P) + 

2K{l + C){Ep + P)-M' 



1 - f M3 

+ 



arctan 



(1 + if AKPT [ \ TM 

2K{l + ^){Ep-P)-M^ 



— arctan 



TM 



.(22) 



In fig. 121 we plot the cross-section as a function of the UHEi/ energy K, for 
neutrino masses m,y ranging from 10"^ eV to 10~^ eV. Each curve corresponds to 
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K [gV] K [gV] 



Figure 2: Cross-section <t^i^{P,K), in cm^, as given by eq. H22() . as a function of the energy 
of the incident neutrino, K, and for neutrino masses of 10~^, 10~^, 10~^ and 10~^ eV (from 
left to right and top to bottom). Each plot shows four curves corresponding to different relic 
neutrino momenta; the dotted (red) curve corresponds to the approximation of relic neutrinos 
at rest, P = 0; the continuous (blue) curve is for Prms = \/{P^ — 6.08 x 10^'' eV and the 
dashed and dot-dashed (black) curves are respectively for 2Pi.nis and 5Prms- 



a particular value of the relic neutrino momentum P. As expected, the thermal 
effects become more and more important as the ratio between the neutrino mass 
and the Cz/B temperature (which we take here as T^q ~ 1.69 x lO"'^ eV) decreases. 
We also compare with the value of the cross-section in the approximation of relic 
neutrinos at rest, (7,^^(0,^), which is obtained by taking the limit of auu{P,K) 
for P ^ 0. As long as the neutrino mass is sufficiently large, that is m,y > 0.01 
eV, the cross-section does not vary significantly over the relevant range of P, 
i.e., the range of momenta selected by the Cz/B distribution function around the 
maximum of P'^fy{P) at Pmax — 3.75 x 10~^ eV. In this case, the position of the 
peak corresponds to the resonant energy K^^^ = M'^/2m^ for target neutrinos at 
rest. For smaller masses, the peak in the cross-section gets broader and shifts to 
smaller UHEz/ energies for increasing relic neutrino momentum. The maximum 
value is approximately constant over the relevant range of momenta, and is about 
one order of magnitude lower than a^uiO, -ft'res)- 

A simpler expression for a^p can be obtained using the mean value theorem by 
evaluating the integrand in eq. ()19p at the midpoint of the integration interval, 
s ~ 2KEp. Taking into account that s+ — s_ = 4KP we get the following 
expression: 
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where all the dependence on K and P is implicit in s = 2K\/m? + P"^. 

We can use these expressions to evaluate the integral in eq. (fTHj) numerically. 
A common approximation consists in neglecting the thermal motion of the back- 
ground neutrinos (see for example jElE!)- The corresponding expression can be 
recovered in our formalism by evaluating the cross-section (jSBI) in P = 0, i.e., 
in So = 2myK, or equivalently by taking the limit for P — of the general 
expression (j221)- The damping reads in this case: 



7°(ir) = a,p(so)n, (24) 

2 Km 

= 2V2GfTM — — — n^, (25) 

with 

poo p 

n, = n,= / —P^f^{P). (26) 
Jo 27r^ 

Substituting the full expression for the total Z width to fermions, eq. (|TT|) 
into eq. (j2SI), one can easily recover the result of [Hj. The expressions of jTHj can 
be obtained by further evaluating the cross-section at the pole of the resonance 
2mKres = M"^ (narrow- width approximation). 

The corresponding results are plotted in fig. El where we directly compare 
the approximations for the damping factor, eq. ()25j) and eq. (jSHI), with the exact 
expression obtained from eqs. ()18j) and ()22j) . The comparison is done for different 
values of the neutrino mass. One notices that the net effect of thermal broadening 
is a reduction of the damping, which affects the transmission probability and the 
depth and shape of the absorption dips. 

Two effects combine: the modification of the cross-section peak due to its 
dependence on Ep, and the presence of the thermal distribution which selects a 
range of relic neutrino momenta comparable to the temperature of the CuB. Both 
effects contribute to smearing out the resonance peak and shifting its maximum 
value from K^es to lower energies. As long as the ratio rriy/Ty is > 100, the 
position of the peak is not significantly affected and the net effect is to reduce the 
damping. For smaller masses, the effect of thermal broadening is much stronger 
and the damping gets spread over a larger range of UHEz/ energies, resulting in 
a worse definition of the absorption dip. 

The comparison of the exact and approximate results for the damping in fig. El 
illustrates that none of the approximations is valid over the range of neutrino 
masses we consider. In the subsequent analysis we therefore keep on working 
with the full expression for the cross-section a^y^K, P), eq. (j221)- 
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Figure 3: Damping factor due io v — v annihilation, in eV, as a function of UHE neutrino 
energy K\ each plot corresponds to a given value of the neutrino mass (10~^ eV, 10~^ eV, 10"'^ 
eV, 10~^ eV from left to right and from top to bottom) and displays the integrated damping, 
7^p, from the full cross-section eq. (black, continued curve) and from the approximated 
cross-section eq. (1^ (blue, dashed curve), as well as the approximation for neutrinos at rest, 
7°p, from eq. (red, dotted curve). The three curves are superposed for rrii, = 0.1 eV, and 
so are the two approximations for rriv — 0.01 eV. 



2.3 Transmission Probability across the relic neutrino back- 
ground 

High energy neutrinos can travel cosmological distances almost without interact- 
ing. To calculate the damping by the Cz/B one has to take the expansion of the 
Universe into account. The standard procedure is to calculate the transmission 
probability by integrating the damping factor over the UHEi/ path, in terms of 
the redshift z, back to the source position z^,. The formula (0) is then generalised 
as follows: 



Pt{Kq,Zs) = exp 



-j,,{Ko{l + z)) 



(27) 



H{z)il + z) 

with the neutrino temperature T^, = Tj,o(l + z). Both Kq and T^o refer to 
the quantities in today's Universe. For the Hubble factor we take H{z) = 
i/o\/0-3(l + + 0.7 as suggested by recent observations |12^, with the numeri- 
cal value of Hq from 23 ; • 

Eq. (j2Zj) encompasses two effects due to the expansion of the Universe. First, 
the UHEz/ energy gets shifted, K Kq{1 + z). This broadens the absorption dip 
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Xo[cV] KolcV] 

Figure 4: Transmission probability Pt{Ko function of the incident neutrino energy as 

detected on Earth, Kq, for an UHE neutrino source located at redshifts Zg = 1, 5, 10, 20 (from 
top to bottom in each plot) and for a neutrino mass rui, = 10~^, 10~^, 10"'^, lO""* eV. The 
continued, black curves corresponds to the full 7i/p from eqs. (fl^ and (|^ . while the dotted 
(red) curves are for the approximation of relic neutrinos at rest, 7^^^, from eq. H25|) . 



even in the approximation of relic neutrinos at rest, since the resonance energy 
changes along the UHEz/ path. Second, the temperature of the Cz/B gets shifted, 
Ti, T^q{1+z), i.e., the thermal bath of relic neutrinos is hotter at earlier times. 
This directly affects the ratio between and T^y and results in a modification 
of the absorption properties of the Cz/B with respect to the UHEz/s that cross it. 

In fig. m we show the transmission probability for an UHEz/ emitted at a fixed 
redshift Zs = 1, 5, 10 or 20, as a function of its present energy, Kq, and compare 
it with the results obtained in the approximation of relic neutrinos at rest. 

As long as my/T^ > 10^, the shape of the absorption dip is not affected by 
thermal broadening and is rather sharply delimited, at high energies, by the bare 
resonant energy for the propagating neutrino, i^'res = M^/(2m,y), and at low 
energies by the redshifted resonant energy i^res/ll + z). Evaluating the position 
of these points would in principle allow us to determine the value of the neutrino 
mass as well as the redshift at which the UHEz/ was emitted. As m^/Tu decreases, 
however, the absorption dips gets shallower and broader, which complicates the 
extraction of and z^. The position of the minimum transmission probability 
is also shifted to lower energies. The effect of thermal motion increases with the 
redshift z since UHE neutrinos from more distant sources are emitted in hotter 
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backgrounds. 



3 Applications 

The results presented so far deal with a monoenergetic source of UHE neutrinos 
located at a given redshift; let us now illustrate the potential effects of thermal 
motion on the process of absorption of UHE neutrinos in two more realistic 
contexts of physical relevance. 

3.1 Absorption lines in a realistic UHE neutrino flux 

The flux of UHE neutrinos at Earth depends on one hand on the mechanisms 
at work in the sources, which determine the injection spectrum of the UHE 
neutrinos, and on the other hand on the spatial and temporal distribution of the 
sources themselves. We follow here the approach of ^HI and express the UHEz/ 
flux in function of the present neutrino energy Kq as: 

1 f°° 7 

J^,{Ko) = — — - Pt{Ko, z) C,{Ko, z), (28) 
Att Jo H{z) 

where C^{Kq,z) is the neutrino source emissivity distribution, given in terms of 
the redshift z and the present energy of the UHEz/. In the hypothesis of identical 
injection spectra for all sources, one can factorize the dependence in the redshift 
and write 

C,iKo,z)=r]{z) J,{Ko), (29) 

where ri{z) describes the distribution of the sources in the Universe, and J^^Kq) 
gives the number of neutrinos emitted per unit of energy by each of these sources. 
We use the following, standard ansatz [ini EEl 121] : 

7]{z) = r/o (1 + zT - z,^^) ^(z^ax - z)- (30) 
UK) = J, K-'^ e{K^^,-K). (31) 

Eq. (j30|l is suitable for an approximate description of UHE neutrino sources 
distribution in models ranging from astrophysical acceleration sites ("bottom- 
up" mechanisms, for which we can take n ~ 4 and Zmax < 10) to exotic, non- 
accelerator sources (which have n ~ 1 to 2 and may extend to a much larger 
2;max)- In both cases, we take the lower bound for the source distribution to be 
Zmin = 0. The spectral index a typically ranges between 1 and 2, depending on 
the production mechanism considered; we assume here that it is constant over the 
range of energies examined and do not consider the possibility of broken power- 
law spectra above the GZK energy. We also suppose that -ft'max > -^res(l + 
in all our analysis. Under these assumptions, as pointed out in jJE], the only 
dependence on the spectral indexes a and n enters through a difference n — 
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a. We consider here, for the purpose of illustrating our results, two distinct 
situations: n — a = 2 which could describe the UHE neutrino flux expected from 
an astrophysical, bottom-up-type source, and — a = which would rather be 
associated to UHE neutrino fluxes produced in top-down processes. Results are 
presented in fig. El which displays the UHEz/ flux, eq. (j2HI), as a function of the 
present energy Kq of the UHEz/, after normalization to the flux in the absence 
of absorption effects, jF^o (obtained by replacing P-y{Kq,z) = 1 in eq. (j2Hl))- 
The first column corresponds to n — a = 2 and for each value of my we show 
several curves corresponding to different redshift limits for the source population, 
-^max = 2, 5, 10, while the second column corresponds to n — a = and redshift 
limits 2;max = 10, 20. 

In that range of parameters, one can see that thermal effects do not signifi- 
cantly affect the shape of the absorption dip in the UHE neutrino flux as long 
as iriy/Ty > 10^. In particular, the endpoint of the dip at high energies, corre- 
sponding to -ft"res, is well-defined and can be used to estimate the absolute value 
of TJiy, provided the absorption dip is not too shallow nor too narrow to be re- 
solved experimentally. The global shape of the dip clearly depends on the value 
oi n — a. Assuming that the the spectral index a can be determined from the 
measurements of the UHEi/ spectrum in the range of energies which is not af- 
fected by the absorption, one could then obtain information on n and therefore on 
the development of the source population. From the endpoint of the absorption 
dip at low energies, which corresponds to the resonance energy of the neutrinos 
emitted at the largest redshift, Kq{1 + Zmax), one can also estimate the epoch at 
which the UHEz/ sources appeared. 

For smaller values of the ratio rriy/Ty, the situation is significantly compli- 
cated by the thermal motion of relic neutrinos. As a result of the broadening of 
the transmission probability, the dips get shallower and can extend over several 
orders of magnitude in energy, depending on the maximum redshift chosen for 
the source distribution. This might complicate their observation, especially at 
small redshifts. From the figure one sees indeed that for z^^^^ < 2, the flux won't 
be depleted more than a 5 %. Results will also be more difficult to interpret in 
terms of a prediction for the neutrino mass and the maximum redshift for the 
population of sources, since the position of the end points of the absorption dip 
is not so clearly defined anymore. On the other hand, maximum absorption is 
now achieved at lower energies Kq. The shift is significant, even more than an 
order of magnitude for very small masses, rriy ^ 10~^ eV. In view of the inverse- 
power-law form of the UHEz/ energy spectrum, this could significantly help to 
improve the detection potential, even though the energy range 10^^ — 10^^ eV, is 
currently beyond the reach of the majority of UHEz/ experiments planned (see 
for example ^16ij). 
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Figure 5: UHE neutrinos fluxes in presence of damping, J^, normalized to the corresponding 
flux in absence of interactions, J-,yo, for a neutrino mass rrii, = 10~^, 10~^, 10~^, 10~^ eV (from 
top to bottom). The black (continued) curves are for the exact expression using eq. (|22(l while 
the dotted, red curves are for the approximation of neutrinos at rest, eq. H25(l . The left column 
is for n — a = 2 and Zmax = 2, 5, 10 (from top to bottom in each plot), while the right column 
is for n — a = and Zmax = 10, 20. 



3.2 Absorption lines due to relic neutrino clustering 

The possibility that massive rehc neutrinos cluster onto dark matter halos has 
been intensively studied, in particular in relation with the generation of the UHE 
cosmic rays through the Z-burst mechanism ^13 IHl E] ■ Recent works have pre- 
sented revised estimations of the density profiles and typical spatial extension of 
the neutrino clusters ||28j|29j, giving overdensities of the order of 10 — 10^ with re- 
spect to the present background density. They can extend on scales L ~ 0.01 — 1 
Mpc, depending on the neutrino mass, the mass of the attracting halo, and its 
velocity dispersion (which is typically of the order of 200 km/s for a galaxy and 
1000 km/s for a galaxy cluster). This last parameter also constrains the epoch at 
which clustering can start, since neutrinos cannot be efficiently trapped as long 
as their mean velocity, 

eV 

~ 1.6 10^ (1 + z) (— ) km s-\ (32) 
m 

is larger than the velocity dispersion of the attracting galaxy or galaxy cluster j29j . 
For neutrinos with masses < 1 eV, clustering will thus take place at very small 
redshifts and we can safely ignore the effect of the expansion of the Universe in 
this analysis. Other important limiting factors to the clustering of neutrinos on 
large scales are Pauli blocking and the limit on the maximum phase-space density, 
as described in j3n|. They actually imply that only neutrinos with mass > 1 
eV will efficiently cluster on galactic halos, on typical scales Lq ~ 50 kpc, while 
neutrinos with mass rrii, > 0.1 eV can cluster on scales as big as Lc ~ 1 Mpc in 
halos associated to (super-) clusters of galaxies [HIl EHl EH] ■ 

The effect of neutrino clustering is limited to small scales and the overden- 
sities are not large enough to have a significant incidence on the damping of 
UHEi/ travelling on cosmological distances, except maybe in the case of galaxy 
superclusters like Virgo |,31j . 

To compute the absorption by clustered neutrinos, one has to substitute a 
suitable distribution function f^^{P) inside the cluster into eq. ^| To make a 
simple estimation of the effect, we assume that f^^{P) does not depend on the 
position, i.e., we take a constant neutrino overdensity. Outside of the cluster, the 
neutrino density is that of the Cz/B. For f^^{P), we make the ansatz of a modified 
Fermi-Dirac distribution 

,1 1 e"*/^" + 1 

= 2 eiP'^y^'^ + l ^^^^ 

This distribution parametrises reasonably well the distribution functions pre- 
sented by j2ni in terms of a single parameter $, keeping the temperature Ti, 
unchanged from the Cz^B. The neutrino density which corresponds to eq. IHHlis 

< = -^(l + e-*/^^)Li3(-e*/^^), (34) 
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Figure 6: Transmission probability for a cluster of extension 1 Mpc, made of neutrinos of 
mass 0.1 eV, with a constant neutrino density n^' — 10'^ no,y and n'^ = 10^ no,^. 




where Li3(x) is the trilogarithm function. For a given overdensity factor N'^^ we 
solve n'^ = N'^^HyQ numerically for $. 

The clusters we consider in the following have an overdensity factor A^'^' be- 
tween 10 and 10^ according to the results mentioned earlier, and a spatial exten- 
sion L^'. We study the case where the cluster is located between the source of the 
UHEz/ and the observer and compute the transmission probability for neutrino 
masses of 1 eV and 0.1 eV. 

As expected, the effect of the thermal motion of the relic neutrinos is in general 
negligible or small due to the relatively small overdensities achievable. We have 
to saturate the bounds on the parameters to get a significant effect, as shown 
in fig. ini which displays the transmission probability with and without thermal 
effects, for a cluster of relic neutrinos with mass = 0.1 eV and extension 1 
Mpc. For a maximal overdensity factor A^'^' = 10^ the thermal motion reduces 
the maximum absorption probability across the cluster from k, 55% to ~ 35%, 
contributing to reducing its effect respect to the non-clustered Cz/B absorption 
probability shown in fig. |31 



4 Conclusions 

Using the formalism of finite-temperature field theory we have calculated the 
damping factor of an UHEz/ propagating through the Cz^B including the effects of 
the thermal motion of the relic neutrinos in a systematic way. This allowed us to 
generalize the expressions for the transmission probability Pt that are commonly 
used in the literature. 

From the exploration of the parameter space allowed by cosmological and as- 
trophysical constraints as well as by current limits on the neutrino mass, we see 
that thermal effects significantly affect the shape and position of the absorption 
dips in a realistic UHEz/ flux as soon as the ratio between the neutrino mass and 
the Cz/B temperature goes below ~ 10^, i.e., well before the relic neutrinos be- 
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come relativistic. As expected, the effect essentially consists in smearing out the 
dip and shifting its minimum value to lower energies. This will complicate the 
observation of the dip in a real experiment measuring neutrino fluxes, especially if 
the UHEi/ source population is concentrated at small redshifts, producing rather 
shallow and extended dips. The shift of the absorption dip to lower energies, 
where neutrino fluxes are expected to be higher, increases in principle the poten- 
tial of detection of the effect respect to the case of a neutrino at rest with the same 
mass. Still, the situation could become more intricate if the pattern of neutrino 
mass eigenstates is such that their combined effect results in a superposition of 
dips with different depths and extensions. 

In the context of neutrino mass spectroscopy, we see from the examples that 
thermal effects do not affect the determination of the endpoint of the absorption 
dip, and hence the possibility of extracting information on the absolute neutrino 
mass, as long as m^, > 0.01 eV, which is verified by at least one neutrino in 
the currently favoured mass schemes [221 • As the absorption dips get broader 
and shallower, the prospects for determining efficiently the resonance energy get 
worse as the endpoint of the absorption dip is no longer sharply defined. 

Finally, as another application of our formalism, we have also investigated the 
transmission probability for UHE neutrinos propagating through a relic neutrino 
cluster. In the standard context of neutrino clustering around galaxies or galaxy 
clusters, we found that the thermal motion can further reduce the absorption 
effect of the cluster. This absorption is small compared to the absorption of 
UHEz/ in the CuB for neutrinos travelling cosmological distances, since clustering 
only occurs at small redshifts. 

Acknowledgments 

We would like to acknowledge support by CONACyT under grants 34868-E and 
46999-F and by DGAPA-UNAM under grants PAPIIT IN116503, IN119405 and 
IN112105. 

References 

[1] p. J. E. Peebles, Principles of Physical Cosmology, Princeton University 
Press, 1993. 

[2] E.W. Kolb and M.S. Turner, The Early Universe, Addison- Wesley, Redwood 
City, 1990. 

[3] T. J. Weiler, Phys. Rev. Lett. 49, 234 (1982). 
[4] T. J. Weiler, Astrophys. J. 285, 495 (1984). 



16 



[5] P. G ondolo, G. Gelmi ni, and S. Sarkar, Nucl. Phys. B 392, 111 (1993) 
|arXiv:hep-ph/9209236|. 

E. Roulet, Phys. Rev. D 47, 5247 (1993). 

S. Yoshida, Astropart. Phys. 2 (1994) 187. 

S. Yoshida, H. y. Dai, C. C. H. Jui, and P. Sommers, Astrophys. J. 479, 547 



[6] 
[7] 
[8 

[9 

[10 

[11 
[12 

[13 

[14 

[15 



(1997) |arXiv:astro-p h/9608186]. 

D . Fargion, B. Mele, and A. Sahs, Astrophys. J. 517, 725 (1999) 
|arXi v:astro-ph/ 9710029|. 

T. J. Weiler, Astropart. Phys. 11, 303 (1999) [arXiv:hep-ph/9710431J . 

E. Waxman, |arXiv:astro-ph/9804023| 

S. Yoshida, G. Sigl, a nd S. J. Lee, Phys. Rev. Lett. 81 (1998) 5505 
[arX iv:hep-p h/9808324l. 

J. J. B lanco-Pillado, R. A. Vazquez, and E. Zas, Phys. Rev. D 61, 123003 
(2000) |arXiv:astro- ph/9902266|. 

H. Pas and T. J Weiler, Phys. Rev. D 63 (2001) 113015 
I arXiv: hep-ph/0101091| . 

D. Fargion, M. Grossi, P. G. De Sanctis Lucentini, C. Di Troia, and 
R. V. Konophch, arXiv:astro-ph/0102426 (in Dark Matter in Astro and 
Particle Physics, Dark 2000, Heidelberg, 10-14 July, Ed. H.V.Klapdor- 
Kleingrothaus, Springer 2001, p. 455-468). D. Fargion, P. G. De Sanctis Lu- 
centini, M. Grossi, M. De Santis, and B. Mele, Mem. Soc. Ast. It. 73 (2002) 
848 | arXiv:hep-ph/0112014, . 

[16] B. Eberle, A. Ringwald, L. Song, and T. J. Weiler, Phys. Rev. D 70, 023007 
(2004) |arXiv:hep -ph/0 401203) . 

[17] G. Barenboim, O. Mena Requejo, and C. Quigg, Phys. Rev. D 71 (2005) 
083002 ar Xivrhep-ph/0412122i . 

[18] J. C. D'Oli vo and J. F. Nieves, Phys. Rev. D 52 (1995) 2987 
I arXi v:hep-p h/9309225 1 . 

[19] J. C . D'O livo and J. F. Nieves, Phys. Lett. B 359 (1995) 148 
|arXiv:hep-ph/9508394|. 

[20] F. Berends et ai, "The Z line shape", in CERN Yellow Report 89-08 (1989) 
Vol. 1, eds. G. Altarelli, R: Kleiss & C. Verzegnassi, p. 89. 



17 



[21] D. Y. Bardin, A. Leike, T. Riemann, and M. Sachwitz, Phys. Lett. B 206, 
539 (1988). 

[22] A. Borrelli, M. Consoli, L. Maiani, and R. Sisto, Nucl. Phys. B 333, 357 
(1990). 

[23] S. Eidelman et al, Phys. Lett. B 592, 1 (2004). 

[24] P. Langacker, J. P. Leveille, and J. Sheiman, Phys. Rev. D 27, 1228 (1983). 

[25] D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 148 (2003) 
175 |arXiv:a stro-ph /0302209 . 

[26] O. E. Kalashev, V. A. Kuzmin, D. V. Semikoz, and G. Sigl, Phys. Rev. D 
66 (2002) 063004 |arXiv:hep-ph/0205050| 

[27] D. V. Semikoz and G. Sigl, JCAP 0404 (2004) 003 |arXiv:hep-ph/0309328| . 

[28] S. Singh and C. P. Ma, Phys. Rev. D 67, 023506 (2003) 
I arXiv:astro-ph/0208419j . 

[29] A. Ringwald and Y . Y. Y. Wong, JCAP 0412, 005 (2004) 
^arXiv:hep-ph/0408241 . 



[30] S. Tremaine and J. E. Gunn, Phys. Rev. Lett. 42 (1979) 407. 

[31] A. Ringwald, T. J. Weiler, and Y. Y. Y. Wong, |a^Xw:astro-ph/ 0505563 



[32] J. F. B eacom and N. F. Bell, Phys. Rev. D 65 (2002) 113009 
farX iv:hep-ph/0204111j . 



18 



